Astronomy & Astrophysics manuscript no. 13234 


©ESO 2010 


August 16, 2010 





The ultracool-field dwarf luminosity-function and space density 
from the Canada-France Brown Dwarf Survey * 

C. Reyle 1 , P. Delorme 2 , C.J. Willott 3 , L. Albert 4 , X. Delfosse 5 , T. Forveille 5 , E. Artigau 6 , L. Malo 6 , G.J. Hill 7 , and R. 

Doyon 6 

1 Observatoire de Besancon, Universite de Franche-Comte, Institut Utinam, UMR CNRS 6213, BP 1615, 25010 Besancon Cedex, 
France 

e-mail: celine@obs-besancon.fr 

2 School of Physics and Astronomy, University of St Andrews, North Haugh, St Andrews KY16 9SS, United Kingdom 

3 Herzberg Institute of Astrophysics, National Research Council, 5071 West Saanich Rd, Victoria, BC V9E 2E7, Canada. 

4 Canada-France-Hawaii Telescope Corporation, 65-1238 Mamalahoa Highway, Kamuela, HI96743, USA 

5 Laboratoire d'Astrophysique de Grenoble,Universite J.Fourier, CNRS, UMR5571, Grenoble, France 
(3JT)' 6 Departement de physique and Observatoire du Mont Megantic, Universite de Montreal, CP. 6128, Succursale Centre- Ville, 

Montreal, QC H3C 3J7, Canada 
7 McDonald Observatory, University of Texas at Austin, 1 University Station C1402, Austin, TX 78712-0259, USA 



o 

(N 



< 



ABSTRACT 



Context. Thanks to recent and ongoing large scale surveys, hundreds of brown dwarfs have been discovered in the last decade. The 
j"y^ Canada-France Brown Dwarf Survey is a wide-field survey for cool brown dwarfs conducted with the MegaCam camera on the 

Canada-France-Hawaii Telescope telescope. 

Aims. Our objectives are to find ultracool brown dwarfs and to constrain the field brown-dwarf luminosity function and the mass 
, function from a large and homogeneous sample of L and T dwarfs. 

£1^ Methods. We identify candidates in CFHT/MegaCam i' and z' images and follow them up with pointed near infrared (NIR) imaging 

on several telescopes. Halfway through our survey we found ~50 T dwarfs and ~170 L or ultra cool M dwarfs drawn from a larger 
sample of 1400 candidates with typical ultracool dwarfs i' - z' colours, found in 780 square degrees. 

Results. We have currently completed the NIR follow-up on a large part of the survey for all candidates from mid-L dwarfs down 



o 

to the latest T dwarfs known with utracool dwarfs' colours. This allows us to draw on a complete and well defined sample of 102 
^ i ' ultracool dwarfs to investigate the luminosity function and space density of field dwarfs. 

Conclusions. We found the density of late L5 to TO dwarfs to be 2.0+%* x 10~ 3 objects pc~ 3 , the density of TO. 5 to T5.5 dwarfs to be 
x 10 3 objects pc~ 3 , and the density of T6 to T8 dwarfs to be 5.3*22 x 10 3 objects pc~ 3 . We found that these results agree better 
with a flat substellar mass function. Three latest dwarfs at the boundary between T and Y dwarfs give the high density 8.3*2° x 10~ 3 
objects pc~ 3 . Although the uncertainties are very large this suggests that many brown dwarfs should be found in this late spectral type 
. range, as expected from the cooling of brown dwarfs, whatever their mass, down to very low temperature. 
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(3 1. Introduction 

y—i ■ * Based on observations obtained with MegaPrime/MegaCam, a joint . 

! • ■ project of CFHT and CEA/DAPNIA, at the Canada-France-Hawaii Brown dwarfs ver y low-luminosity objects. Even at their flux 

Telescope (CFHT), which is operated by the National Research Council maximum in the near and mid-infrared, brown dwarfs are more 

(NRC) of Canada, the Institut National des Sciences de l'Univers of the than ten magnitudes fainter than solar-type stars. That explains 

Centre National de la Recherche Scientifique (CNRS) of France, and the very low number of detected substellar objects compared to 

■ the University of Hawaii. This work is based in part on data products that of known stars, although they probably represent a sizeable 
produced at TERAPIX and the Canadian Astronomy Data Centre as fraction of the stellar population in our Galaxy, 
part of the Canada-France-Hawaii Telescope Legacy Survey, a collab- The first discove ries of brown dwarfs were made by 

T^fiSK*?*? ™t ^ NRS ' Ba l ed T ° n *™ i0nS made T3 wi * Naka iima et all (fl995h aro und the nearb y star Gliese 229 and 

the ESO New Technology Telescope at the La Silla Observatory. Based , ITT^'Tr-^'T'iTTTTTrn, ilTr - ; — 1 : — n7T7UT?i . , . , 

, .. , t . j . i, . . ™ . ,. , . by Staufferet al. (1994) and Rebolo et al. 1995) in the Pleiades. 

on observations obtained at the Gemini Observatory, which is oper- i ' - — -n' ' j 'r- — . 1 ,1 TTZZZt — , | ■' , . , n 

ated by the Association of Universities for Research in Astronomy, Delfosse_etalJ dl997|); | Ruiz et al.| dl997|) and | Kirkpatnck et al. | 

Inc., under a cooperative agreement with the NSF on behalf of the U222II found the first field brown dwarfs. Since then, several 

Gemini partnership: the National Science Foundation (United States), hundreds of field brown dwarfs have been identified, most of 

the Science and Technology Facilities Council (United Kingdom), 

the National Research Council (Canada), CONICYT (Chile), the Foundation. Based on observations made with the Nordic Optical 

Australian Research Council (Australia), CNPq (Brazil) and CONICET Telescope, operated on the island of La Palma jointly by Denmark, 

(Argentina). Based on observations with the Kitt Peak National Finland, Iceland, Norway, and Sweden, in the Spanish Observatorio del 

Observatory, National Optical Astronomy Observatory, which is op- Roque de los Muchachos of the Instituto de Astrofisica de Canarias. 

erated by the Association of Universities for Research in Astronomy, Based on observations made at The McDonald Observatory of the 

Inc. (AURA) under cooperative agreement with the National Science University of Texas at Austin. 
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them thanks to large-scale optical and near-infrared imaging 
surveys because they can be identified by their red optical mi- 
nus near-infrared colours. Most of them have been found in 
the Two-Micron All-Sky Survey (2MASS, Skrut skie et ajl2 006) 
and in the Sloan Digit al Sky Survey (SPSS. lYork et al.l l2000). 
See[Kirkpatrick] d2005l) for a review on brown dwarf discoveries. 

The new generation of large-area surveys uses deeper im- 
ages. As a consequence, the number of known brown dwarfs 
increases and new types of rarer and fainter brown dwa rfs are 
detected (iDelorme et al.ll2008at iBurnin gham et al. 2008). Such 
surveys are the UKIR T Infrared Deep Sky Survey (UKIDSS, 
lLawrence et al.l 120071) and the one we undertook, the Canada - 
France-Brown-Dwarf Survey (CFBDS, IDelorme et a l. 2008b). 
The Canada-France-Brown-Dwarf Survey is based on deep 
multi-colour MegaCam optical imaging obtained at the Canada- 
France-Hawaii Telescope (CFHT). We expect that complete 
characterisation of all our candidates will yield about 100 T 
dwarfs and over 400 L or very late-M dwarfs, which will ap- 
proximately double the number of known brown dwarfs with a 
single, well-characterised survey. 

With that large number of identified brown dwarfs, it be- 
comes possible to define uniform and well-characterised sam- 
ples of substellar objects to investigate their mass and luminos- 
ity functions. The knowledge of these functions is essential in 
several studies. In terms of Galactic studies, it gives clues on the 
baryonic content of the Galaxy and contributes to determine the 
evolution of the Galaxy mass. In terms of stellar physics studies, 
it gives constraints on stellar and substellar formation theories. 

It has been shown that a single log-normal functio n could fit 
the m ass function from field stars to bro wn dwarfs (see lChabrierl 
120031) ). (see also iLuhman et al.l |2007)). Our Galaxy probably 
counts several 10 10 of brown dwarfs, hundreds of them in the 
Solar neighbourhood! 

To validate these assumptions, which are mainly based on 
the study of brown dwarfs in stars clusters, it is necessary to first 
refine the field brown dwarf luminosity-function. An initial esti- 
m ate of the local sp ace density of late T dwarfs have been made 
bvlBurgasserl (l2002R from a sample of 14 T-dwarfs. Allen et al. 
(2005) used these results combined with data compiled from a 
volume-limited sample of late M and L dwarfs (ICruz et al.l2 003) 
to compute a luminosity function. Recently, a detailed investiga- 
tion on a volume-li mited sample of fi eld late-M and L dwarfs has 
been performed bv lCruz et al.l d2007l) . A similar e mpirical inves- 
tigatio n of field T-dwarfs have been presented bv lMetchev et al.l 
(2008). Still the number of objects in their sample is relatively 
small (46 L-dwarfs and 15 T-dwarfs, respectively) and the field 
brown dwarf luminosity-function remains poorly constrained. 
Thus further efforts are needed to measure the space density of 
brown dwarfs. 

At mid-course of the CFBDS survey, we are able to de- 
fine an homogeneous sample of 102 brown dwarfs redder than 
i' - z' =2, from the mid-L dwarfs to the far end of the brown 
dwarfs observed sequence at the T/Y transition, and to derive a 
luminosity function. These objects are drawn from a 444 square 
degree area (57% of the total area) where all candidates with 
i' - z' > 2.0 have been followed-up with near-infrared pho- 
tometry and, for the reddest of them, with spectroscopy. They 
represent a sub-sample in colour (/' — z' < 2.0) and in magni- 
tude (z' < 22.5) of the hundreds of brown dwarfs found with 
i' - z' > 1 .7 on the entire survey. Sect. [2] briefly describes 
the CFBDS and the related observations. The construction of 
a complete and clean sample of L5 and cooler dwarfs is ex- 



plained in Sect. 13.21 In Sect. [5] we compare our results with 
prior studies and link them to the mass function of brown dwarfs. 
Sect. [4] presents the luminosity function derived from this sam- 
ple. Conclusions are given in Sect. [6] 

2. Observations 

The brown dwarf sample is d rawn from the CFBDS . The 
survey is fully described in IDelorme et al.l d2008bi) . The 
CFBDS is a su rvey in the i' and z' filters conducted 
with MegaCam dBoulade et al l 120031) at the CFHT. It is 
based on t wo existing surveys, the CFHT Legacy Survey 
(CFHTLSEL ICuillandre & Bertinl 120061) a nd the Red-sequence 
Cluster Survey (RCS-2. lYee et al.l 120071) . complemented with 
significant Principal Investigator data at CFHT. The survey is 
also extremely effective at finding high-redshift quasars. The 
parallel programme is called the C anada-France High-z Quas ar 
Survey, and results are presented in lWillott etaE (2007. 20091). 

The i' and z' imaging part of the 780 deg 2 CFBDS is nearing 
completion to typical limit of z' = 22.5, probing the brown dwarf 
content of the Galaxy as far as 215 pc for the mid-L dwarfs, 
180 pc for the early-type T dwarfs and 50 pc for the late-type T 
dwarfs. 

The reddest sources are then followed-up with pointed J- 
band imaging to distinguish brown dwarfs from other astronom- 
ical sources, and spectra are obtained for the latest type dwarfs. 

Throughout this paper, we use Vega magnitudes for the J- 
band and AB magnitudes for the i' and z' optical bands. 



2. 1 . Optical imaging 

We only briefly describe the processing and photometry of the 
imaging observations in this p aper. A full description can be 
found in IDelorme et al.l d2008bl) . 

Pre-processing of the MegaCam images is carried out at 
CFHT using the ELIXIR pipeline. This removes the instrumen- 
tal effects from the images. We then run our own algorithms 
to improve the astrometry and check the photometry. Finally, 
we stack the images (if there is more than one exposure at a 
given position) and register the images in the four different fil- 
ters. Photometry is carri ed out with an adaptation of SExtractor 
( Bert in & Arno uts 1996), which uses dual-image multiple point- 
spread function fitting to optimise the signal-to-noise ratio (S/N) 
of point sources. 

Candidate brown dwarfs and quasars are initially identified 
on the optical images as objects which have very high i' - z' 
colours. A 10 cr detection in z' is required but no constraint is 
set on i'\ in several cases, the targets are i' -dropout s and only 
i' — z' lower limits can be determined. As shown in lWillott et al.l 
(2005), photometric noise causes many M dwarfs with intrinsic 
colours of 0.5 < i' — z 1 < 1 .5 to be scattered into the region of 
the diagram at i' - z' > 1 .5 where we would expect to find only 
L or T dwarfs and quasars. The huge number of these M stars 
would require a lot of telescope time for complete follow-up. 
Therefore, we limit our survey to objects observed to be redder 
than i' - z' > 1 .7. Changing this criterion reduces the number of 
M dwarf contaminants by ~ 95%. 

2.2. Near-infrared imaging 

As shown by iFan et al] (1200 ll) . brown dwarfs and quasars 
(and artefacts) can be separated with NIR /-band imaging. 



http://web.mit.edu/ ajb/www/thesis 



http://www.cfht.hawaii.edu/Science/CFHLS/ 
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IWillott et alJ (I2005h described in detail the method for identify- 
ing high-redshift quasars and brown dwarfs using MegaCam op- 
tical plus near-IR imaging. The very red i' - z' of high-redshift 
quasars is caused by deep Lyman-ff absorption on a relatively 
flat intrinsic spectrum, and they appear significantly bluer than 
brown dwarfs in z' - /. On the contrary the spectral distribu- 
tion of brown dwarfs keeps rising into the / band. We there- 
fore carried out NIR imaging at several observatories: La Silla 
(New Technology Telescope, 3.6m), McDonald (2.7m), Kitt 
Peak (2.1m), La Palma (Nordic Optical Telescope, 2.5m). We 
obtained the / magnitude for all candidates, dwarfs or quasars. 
For dwarfs, the signal-to-noise ratio is 10 to 50 cr. 

Besides pinpointing the few high-redshift quasars that 
contain important clues on the reionization of the Universe 
(IWillottet al.ll2005l l2007l l2009h . the /-band photometry very 
effectively rejects any remaining observational artefacts, as 
well as the more numerous M stars scattered into the brown 
dwarf/quasar box by large noise excursions. 

2.3. Spectroscopy 

Follow-up spectroscopy of the T dwarf candidates with the red- 
dest z' - J colour was then carried out at Gemini. Before the 
GNIRS incident at Gemini-South on April 2007 (see Gemini 
Observatory call for proposal archive, semester 2OO7B0), it was 
used in its cross-dispersed mode to obtain a 0.9-2.4 micron cov- 
erage. Then NIRI at Gemini-North was used in a two-step ap- 
proach: 1) //-band spectra were obtained to determine spectral 
types and 2) T6 and cooler objects were targeted for additional 
coverage in the / band and occasionally in the K band. 

To date we obtained spectroscopic follow-up for 40 T- 
d warfs. Four more T- d warf candidates h a ve been followed-u p 
bv lKnanp et all (l2004h : lchiu et all (120061 ): I Warren et al.1 (l2007h . 
Spectral types range between TO and the T/Y transition with 
6 dwarfs being of type T7 or later. Four of these have already 
been presented in dDelorme et al.ll2008allbl) . Spectral indices of 
the others are given in Table Q] and ar e fully described in a forth- 
coming paper (see also lAlbert et al.ll2009l) . For all of them we 
use an absolute magnitude versus spectral type relation instead 
of the absolute magnitude versus colour used for the other ob- 
jects, (see Sect. 14. 11 1. 



We used public ly available spectra from the L and T 
dwarf data archivfl (iMartfn et al.l | l999tlKirkpatrick et alfc oOO: 
i Geballe et ap 1200 it iLeggett et al.ll2002t iBurgasser et al.li 2003 ; 
Knapp et aTT l2004 iGolimowski et al.l l2004: Chi u et al.ll2006h of 
45 brown dwarfs with spectral types LI to T8. Spectral types for 
L and T dwarfs are gi ven according t o the infrar ed classification 
schem e described in Geball e~et al.1 (120021) and feurgasse r et al.1 
(2006), respectively. 

Thus synthetic colours in the M egaCam filters are computed 
in the AB system (Fukugita et al. 1996) using detector quantum 
efficiency and transmission curves for the atmosphere, telescope, 
camera optics, and filters. We similarly synthesised /-band pho- 
tometry for each of the instruments and / filters used in the /- 
band follow up. These instruments have significantly different 
response curves, which must be taken into account to obtain ho- 
mogeneous selection criteria. 

Synthetic colours are given in Table [2] Figure Q] displays the 
resulting colour-colour diagram. The synthetic colours of each 
brown dwarf are represented by filled squares. The dashed line 
shows the resulting colour-colour relation. Spectral types are 
given along this track, indicating the averaged colour of brown 
dwarfs at a given spectral type. 
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3. Defining a homogeneous, complete, and clean 
sample 

3.1. Photometric classification 

One square-degree MegaCam image contains several hundred 
thousand objects, of which at most a few are brown dwarfs. We 
thus need to strike a careful balance between sample complete- 
ness and contamination. To tune this compromise we need a pre- 
cise knowledge of the colours of brown dwarfs for the exact in- 
struments and filters used in our survey. 

The spectral energy distribution obtained from several pho- 
tometric bands is characteristic of the spectrum of one object. In 
order to get information on the spectral type of the CFBDS can- 
didates based on its photometry, we determined the colours in 
the photometric system used by CFBDS of brown dwarfs with 
known spectral type. 



3 http://www.gemini.edu/sciops/observing-with-gemini/previous- 
semesters/cfp- archive ?q=node/l 1034 



Fig. 1. z' — J synthetic colour versus V — z' synthetic colour 
computed from available spectra in the literature, i' and z' are 
computed in the CFHT/MegaCam filters, / is computed for the 
NTT/SOFI J s filter. The filled squares show the colours of each 
brown dwarf. The dashed line shows the mean colour-colour re- 
lation. Spectral types are given along this track, indicating the 
averaged colour of brown dwarfs with a given spectral type. 
Crosses with error bars show the T dwarfs for which we obtained 
spectroscopic observations. An arrow indicates no detection in 
the /'-band, meaning that the i' - z' colour is a lower limit. 



Owing to the high dispersion of brown dwarf colours, the 
classification of one candidate based on its photometry is not re- 
liable. However, the colours can be used to classify the object 
within a large category such as early L or late L, the overall clas- 
sification of the objects being statistically significant. Table [3] 



4 http ://w w w.j ach. hawaii . edu/~ skl/LTdata. html 
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Table 1. Spectral indices for the the unified scheme of lBurgasser et al.l (l2006|) and spectroph otometric distances. Each index trans- 
lates to a spectral type (between brackets) through linear interpolation using lBurgasser et alJ s Table 5. The adopted spectral type is 
a straight average of the spectral indices that could be measured (upper limits are rejected). 



Designation 


H 2 o-y 


CH 4 - J 


H 2 0-H 


CH 4 - H 


CH 4 - K 


Adopted SpT 


CFBDS J005910-01 1401 


0.037[T8.66] 


0.175[T8.58] 


0.125[>T8] 


0.086[T8.81] 


0.111[>T6] 


T9.0 


CFBDS J025401- 182529 


0.160[T6.40] 


0.350[T6.12] 


0.267[T6.88] 


0.244[T7.06] 


0.163[T6.35] 


T6.5 


CFBDS J025558- 173020 






0.580[T1.29] 


0.980[T0.66] 




T1.0 


CFBDS J025718-124853 






0.412[T4.30] 


0.774[T3.13] 




T3.5 


CFBDS J025805- 145534 


0.555[<T2] 


0.642[T2.32] 


0.568[T1.45] 


1.013[<T0] 


0.723[T1.22] 


T1.5 


CFBDS J025 840- 182648 






0.559[T1.59] 


0.702[T3.49] 




T3.0 


CFBDS J030 130- 104504 






0.325[T5.90] 


0.378[T5.86] 




T6.0 


CFBDS J030225- 144125 






0.306[T6.23] 


0.434[T5.38] 




T5.5 


CFBDS J030226-143719 


0.381[T3.99] 


0.595[T2.83] 


0.408[T4.37] 


0.563[T4.31] 




T4.0 


CFBDS J090139+174051 






0.444[T3.55] 


0.630[T3.85] 




T4.0 


CFBDS J090449+ 165347 






0.408[T4.37] 


0.596[T4.03] 




T4.0 


CFBDS J092250+ 152741 






0.279[T6.68] 


0.238[T7.12] 




T7.0 


CFBDS J102841+565401 


0.043[T8.53] 


0.257[T7.34] 


0.179[T8.35] 


0.147[T8.04] 




T8.0 


CFBDS J104209+580856 


0.160[T6.40] 


0.349[T6.16] 


0.248[T7.21] 


0.269[T6.83] 


0.167[T6.25] 


T6.5 


CFBDS J145044+092108 






0.483[T2.68] 


0.690[T3.55] 




T3.5 


CFBDS J145847+061402 






0.552[T1.69] 


0.871[T2.41] 




T2.0 


CFBDS J145935+085751 






0.405[T4.42] 


0.528[T4.60] 




T4.5 


CFBDS J 150000- 182407 


0.335[T4.89] 


0.441[T5.10] 


0.387[T4.72] 


0.485 [T4.96] 


0.429[T3.61] 


T4.5 


CFBDS J15 1803+07 1645 






0.535[T1.93] 


0.788[T3.06] 




T2.5 


CFBDS J152514+111833 






0.500[T2.43] 


0.833[T2.73] 




T2.5 


CFBDS J152655+034536 


0.464[T2.86] 


0.532[T3.80] 


0.409[T4.35] 


0.574[T4.22] 


0.305[T4.41] 


T4.0 


CFBDS J203737- 192202 






0.642[T0.29] 


0.998[T0.08] 




TO.O 


CFBDS J203841-185012 






0.499[T2.45] 


0.781[T3.10] 




T3.0 


CFBDS J204803-183212 






0.381[T4.82] 


0.537[T4.53] 




T4.5 


CFBDS J2 12243 +042942 






0.571[T1.42] 


0.864[T2.46] 




T2.0 


CFBDS J2 12702+002344 






0.450[T3.32] 


0.732[T3.34] 




T3.5 


CFBDS J214139-033739 






0.529[T2.02] 


0.994[T0.21] 




T1.0 


CFBDS J223856+034947 






0.561[T1.55] 


0.848[T2.60] 




T2.0 



gives the locus of the spectral classes in the colour-colour dia- 
gram. 

In the L-dwarf domain, bright L dwarfs discovered in SDSS 
or 2MASS and with known spectral type that are also identi- 
fied on the CFBDS images allow us to validate the colour-based 
classification. In the T-dwarf domain, we compared the synthetic 
colours as a function of spectral type with the real colours of 
T-dwarf candidates that we followed-up spectroscopically. Their 
locus in the colour-colour diagram is shown by crosses with error 
bars in Fig.Q] While the agreement is good for early T-dwarfs, it 
is poor for the late T-dwarfs where the observed z'—J colours are 
bluer than the synthetic ones. In particular, some of the coldest 
known object are located in the same region of the colour-colour 
diagram as mid-T dwarfs. 

Several reasons can be invoked to explain the discrepancy 
between observed and synthetic colours of late-T. First, the num- 
ber of available spectra of late T-dwarfs is very small (only five 
with spectral type later than T6). Because the colour dispersion 
is quite high for a given spectral type, the mean synthetic colour 
can be biased by one peculiar object. Next, the synthetic colours 
of late T-dwarfs may be uncertain, because for these extreme 
objects with very narrow flux peaks a small error in the opti- 
cal transmission can lead to large colour uncertainties. However, 
as we performed spectroscopic follow-up for all candidates with 
late-T colours, we used their spectral type for classification in- 
stead of colours. 

3.2. The sample 

The CFBDS is composed of several patches - contiguous areas 
on the sky - with sizes ranging from 9 to 79 deg 2 . The galactic 



latitude over a patch is nearly constant, as well as the stellar den- 
sity and the interstellar reddening. The reddening is low for all 
patches (~ 0.01 1 ± 0.009 for RCS2 patches and ~ 0.020 + 0.009 
for CFHTLS patches). Therefore each patch constitutes a sample 
of homogenous data. 

Because the priority was given for the reddest candidates 
(/' - z' — 2 or redder) for /-band follow-up, we can now build 
a complete sample of late L and T dwarf candidates over 16 
patches with a total effective area of 444 deg 2 (Table @). The 
sample contains all candidates with z' < 22.5 and i' — z' > 2.0 
detected in this area. 



This sample contains 249 objects that are candidates cooler 
than L5 on the basis of their i' - z' colour. However, due to the 
photometric errors, many contaminants with true colour i' -z' < 
2.0 enter this sample, such as the more numerous M and early 
L dwarfs. A classification of the candidates is made on the ba- 
sis of their position in the i' - z'/z' - J diagram, as shown in 
Fig. [2] and explained in Sect. 13.1 1 (see Table [3] for a summary). 
Based upon this classification, 102 of the 249 objects remain L5 
and later dwarf candidates (Table 3), the others are dwarfs ear- 
lier than M8, artefacts, or quasars. Dwarfs of spectral type > M8 
to L5 are theoritically not retained by the i' - z' > 2.0 criteria. 
But they could statistically contaminate our >L5 sample if their 
i' - z' colour is reddened by photometric errors; then the z'—J 
criteria will not reject them because they have similar colours 
as the >L5 dwarfs in this wavelength bandpass. This contamina- 
tion has to be carefully estimated (Sect. l3~4l i. The distribution of 
objects within the spectral classes is given in Table|6] 
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Table 2. Synthetic colours computed on publicly available spec- 
tra. 



Name 


spT 


i' 


z' 


J 


i'-z' 


z'-J 


2MASS0345+25 


LI 


18.47 


16.89 


13.84 


1.58 


3.05 


2MASS 1439+19 


LI 


16.93 


15.36 


12.65 


1.57 


2.71 


2MASS0028+15 


L3 


21.51 


19.87 


16.63 


1.64 


3.24 


2MASS2224-01 


L3.5 


19.01 


17.29 


13.88 


1.72 


3.41 


DENIS 1058-15 


L3 


18.87 


17.33 


14.12 


1.54 


3.22 


GD165B 


L3 


20.51 


18.89 


15.63 


1.62 


3.26 


2MASS0036+18 


L4 


17.25 


15.66 


12.31 


1.59 


3.35 


LHS102B 


L4.5 


17.82 


16.04 


12.99 


1.78 


3.05 


SDSS0835+19 


L4.5 


21.12 


19.14 


16.06 


1.98 


3.08 


2MASS1507-16 


L5.5 


17.44 


15.59 


12.69 


1.85 


2.90 


DENIS0205-11 


L5.5 


19.48 


17.36 


14.40 


2.12 


2.96 


SDSS0107+00 


L5.5 


20.98 


18.78 


15.71 


2.20 


3.07 


SDSS0206+22 


L5.5 


20.88 


19.45 


16.42 


1.43 


3.03 


SDSS 1326-00 


L5.5 


21.18 


19.03 


16.15 


2.15 


2.88 


SDSS0539-00 


L5 


18.95 


16.87 


13.85 


2.08 


3.02 


SDSS 1446+00 


L5 


20.72 


18.69 


15.55 


2.03 


3.14 


2MASS0825+21 


L6 


20.39 


18.16 


14.88 


2.23 


3.28 


SDSS0236+00 


L6.5 


21.33 


19.13 


15.99 


2.20 


3.14 


DENIS 1228-15 


L6 


19.44 


17.57 


14.28 


1.87 


3.29 


2MASS 1632+19 


L7.5 


21.25 


18.88 


15.75 


2.37 


3.13 


2MASS2244+20 


L7.5 


21.65 


19.65 


16.32 


2.00 


3.33 


2MASS 1523+30 


L8 


21.53 


19.18 


15.94 


2.35 


3.24 


SDSS0857+57 


L8 


20.28 


17.86 


14.78 


2.42 


3.08 


SDSS1331-01 


L8 


20.59 


18.66 


15.32 


1.93 


3.34 


2MASS0310+16 


L9 


21.34 


18.89 


15.83 


2.45 


3.06 


2MASS0908+50 


L9 


19.85 


17.48 


14.40 


2.37 


3.08 


2MASS0328+23 


L9.5 


21.78 


19.37 


16.34 


2.41 


3.03 


SDSS0805+48 


L9.5 


19.74 


17.96 


14.67 


1.78 


3.29 


SDSS0830+48 


L9 


20.66 


18.26 


15.21 


2.40 


3.05 


SDSS0151+12 


T0.5 


22.49 


19.56 


16.24 


2.93 


3.32 


SDSS0423-04 


TO 


19.73 


17.38 


14.29 


2.35 


3.09 


SDSS0837-00 


Tl 


23.86 


20.43 


16.89 


3.43 


3.54 


SDSS0758+32 


T2 


21.46 


18.49 


14.77 


2.97 


3.72 


SDSS 1254-01 


T2 


21.68 


18.11 


14.65 


3.57 


3.46 


SDSS1214+63 


T3.5 


23.20 


20.24 


16.16 


2.96 


4.08 


OUjjl /JU+1 / 


1 j.j 


11 

ZZJ7 


1 Q 7Q 

17. ly 


1 A 1 1 
10. 1j 


Z.oU 


J. OO 


SDSS1021-03 


T3 


22.12 


19.42 


15.87 


2.70 


3.55 


2MASS0559-14 


T4.5 


21.10 


17.63 


13.56 


3.47 


4.07 


SDSS0926+58 


T4.5 


22.16 


19.25 


15.46 


2.91 


3.79 


SDSS 1624+00 


T6 


22.41 


19.34 


15.19 


3.07 


4.15 


SDSS 1346-00 


T6.5 


22.90 


19.54 


15.48 


3.36 


4.06 


2MASS0727+17 


T7 


22.73 


19.47 


15.18 


3.26 


4.29 


G1570D 


T7.5 


23.20 


19.26 


14.75 


3.94 


4.51 


2MASS0415-09 


T8 


22.97 


19.80 


15.29 


3.17 


4.51 



Table 3. Definition of spectral classes from the i' - z' and z'-J 
colours. 



Spectral class 



• z range 



/ range 



<M8 

M8-L4.5 

L5-T0 

T0.5-T4.5 

>T4.5 



< 2.0 
2.0-2.6 
> 2.6 



<2.5 

> 2.5 

> 2.5 
2.5-3.: 

> 3.8 



3.3. Completeness 

Owing to photometric errors, some objects with true colours 
within our selection criteria are spread out of our selection box 
and are not identified. Given the depth of CFBDS, one can- 
not estimate the completeness with a reference sample of ob- 
jects detected in the images with previously known magnitude 
and colours. To obtain a reference sample we built a point 



Table 4. Patches used to build a complete sample. All candidates 
within our selection box (z' < 22.5 and i' — z' > 2.0) in these 
patches have /-band photometry. 



patch 


survey 


alpha 


delta 


1 


b 


area 






hh:mm) 


(dee) 


(dee) 


(dee) 


(dee 2 


0047 


RCS-2 


00:47 


+00 


120 


-63 


47 


0310 


RCS-2 


03:10 


-16 


201 


-56 


56 


0357 


RCS-2 


03:57 


-07 


197 


-42 


27 


1040 


RCS-2 


10:55 


+57 


149 


+54 


30 


0112 


CFHTLS 


00:45 


+34 


121 


-29 


18 


1237 


CFHTI S 

V J 111 m—JVj 


08:30 


+ 19 


206 


30 


9 


1319 


CFHTLS 


09:15 


+ 17 


213 


39 


15 


1512 


CFHTLS 


10:30 


+ 11 


232 


53 


14 


1514 


RCS-2 


15:14 


+06 


7 


50 


48 


1645 


RCS-2 


16:45 


+40 


64 


41 


28 


2017 


CFHTLS 


13:35 


-08 


231 


53 


13 


2143 


RCS-2 


21:43 


-01 


55 


-38 


79 


2265 


CFHTLS 


15:00 


-17 


342 


36 


14 


3071 


CFHTLS 


19:30 


-22 


17 


-18 


13 


3303 


CFHTLS 


22:10 


-18 


38 


-52 


14 


3400 


CFHTLS 


22:30 


+03 


69 


-44 


19 




late T \ 
► > 



, > early T 

> > 



Fig. 2. z'-J versus i' -z' diagram of the brown dwarf candidates 
in our sample with complete /-band follow-up. Open circles: 
dwarfs earlier than M8. Filled circles: M8 to L dwarfs. Open tri- 
angles: early T dwarfs. Filled triangles: late T dwarfs. The dotted 
line shows the colour-colour relation derived in Sect. 13.11 The 
solid lines show the different spectral type regions defined in 
Table [3] Note that to build the T dwarfs luminosity-function we 
did not use this colour-based classification but used the spectral 
type obtained with spectroscopic follow-up whenever available. 



spread function (PSF) model for each science image using PSfex 
(E. Bertin, private communication) and created fake star s with 
1.2 < i' - z' < 4 and 20 < z' < 23.5 with Skymaker (iBertinl 
120091) . These artificial stars are thus a PSF model scaled to the 
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Table 6. Preliminary classification of the sample based on the 
i' - z' and z' - J colours. The spectral class definitions are given 
in Table [3] Note that we expect most contamination to occur 
in the late L class from the M8 to early L class. The limiting 
distance of detection is also given. 



Object type 


number 


percentage 


distance 


late T dwarfs 


13 


5% 


60 pc 


early T dwarfs 


28 


11% 


100 pc 


late L dwarfs 


61 


25% 


200 pc 


<M8 dwarfs 


77 


31% 


500 pc 


artefacts and quasars candidates 


70 


28% 




total 


249 







desired magnitude, affected with a realistic photon noise. This 
model is oversampled so that it handles correctly the randomised 
non-integer pixel position of the fake stars. We then added these 
stars to the science images from which the PSF model was de- 
rived, which means that any analysis of these fake stars is af- 
fected by the same background noise, bad pixels, saturated stars, 
and cosmic ray distribution as the actual stars are. 

To avoid any significant increase of the image crowding 
caused by this addition, we injected a number of artificial stars 
that is lower than 5% of the number of astrophysical sources on 
the images. This left us with a statistically comfortable sample of 
1 500 000 artificial ultracool dwarfs. Drawing upon this sample 
we were able to derive CFBDS's completeness. 

The resulting images, containing both astrophysical sources 
and fake stars, went through the same analysis and selection 
pipelines used to select true candidates. The injection of fake ob- 
jects in the science images allows us to take into account all ef- 
fects such as cosmic rays, bad pixels, resampling noise, and pos- 
sible selection biases when going through the selection pipeline. 
The effect of faint sources being lost due to the presence of a 
nearby bright star is also perfectly taken into account without 
the use of the somewhat arbitrary masking of a given area around 
saturated stars: their recovery rate directly depends on their mag- 
nitude and separation to the saturated stars actually present in the 
science images. The injection of stars into the actual science im- 
ages is a direct and extremely robust measure of completeness 
regardless of what actually causes the loss of completeness. 

Thus we can count the number of detected objects as a func- 
tion of magnitude on the detection image (z'-band) and colour. 
This is done separately for each patch of the CFBDS. The com- 
pleteness is given by the fraction of recovered fake objects at the 
end of the analysis process. Note that the measured colours and 
magnitudes of fake objects are different from the injected colours 
because of noise and measurement errors. Figure [3] shows the 
completeness averaged over the patches of the RCS-2 survey as 
a function of "true" (that is before going through the analysis 
process) magnitude and colour. We note cy the completeness in 
each bin (/, /') of 0. 1 mag in magnitude and colour. Note that this 
completeness does not merely reflect the number of objects de- 
tected, but the number of detected objects which went through 
our whole pipeline and were selected as ultracool dwarfs can- 
didates with a star-like shape, a signal-to-noise ratio over 10, a 
measured i' - z' colour over our selection criterion (/' - z' > 2.0 
for the sub-sample studied here), and a measured z' magnitude 
below 22.5. 

Moreover, the measured magnitudes are compared with the 
"true" ones to derive an error probability distribution (Fig|4|l. It 
gives the probability for one object in a given patch to have a 
different measured photometry by a given number of <x from its 




20 21 22 

2 MaqnitudefAB) 



Fig. 3. Completeness of the RCS2 component of the CFBDS 
survey as a function of z' magnitude and i' -z' colour. The colour 
bar ranges from 0% (blue, all objects went out of the selection 
limits) to 100% (red, all objects detected). The volume-weighted 
average completeness over our selection criterion (black rectan- 
gle, i' — z' > 2.0 and z' < 22.5) is 73%. This plot also highlights 
the contamination of our sample by objects whose actual colours 
and magnitudes lay outside our selection box but are scattered 
inside by photometric errors. 



true photometry. The result is very close to a Gaussian distribu- 
tion, the only difference is found in the wings at more than 5 <r, 
where the number of deviations stands at 1-2%, several orders of 
magnitude above a Gaussian distribution. This is due to blended 
objects with galaxies, bad pixels, etc. This distribution, which 
also takes into account the correlation between pixels caused 
when resampling any images, will be used when computing the 
contamination of the sample (Sect. l3~4t . 




-4-2 2 



sigma 

Fig. 4. Probability error distribution for one of the RCS2 
patches, in the i' and z' bands. The errors are given in number 
of <t. They are computed by comparing the magnitudes of the 
fake cool dwarfs when injected in the science images with their 
measured magnitudes. 
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zf - J— 23 to zf - 7=2.9, we are confident that possible contami- 
nation from earlier than M8 dwarfs part of our 249-strong brown 
dwarf candidates sample into our zf - J confirmed brown dwarf 
sample trough 4 to 8 <x noise excursion is marginal, (iii) Our 
sample is z' selected, so the flux limit bias or Eddington bias that 
occurs due to a scattering of sources fainter than our zf magni- 
tude limit into the sample will act to increase the number of high 
i' - z! objects, but decrease the number of high zf - /' objects. 
In practice this bias is very small in both directions because we 
have a 10 cr zf limit, (iv) The main reason the scattering from 
low i' - zf to high i' - zf is so bad is because the number den- 
sity of objects with true low if - zf is so much higher than those 
with true high i' —zf. The opposite is actually true for the L and T 
dwarfs: in our i' -zf selected sample, there are more objects with 
true high zf — J (L and T dwarfs) than low zf — J (z~6 quasars). 

The zf — J cut cannot distinguish ultracool dwarfs with true 
colour 1.3 < i' —zf < 2.0 because they have the same zf—J colour 
as late L and early T dwarfs. We therefore need to evaluate the 
number of potential contaminants that populate this colour range 
and can leak into our selection through photometric errors. The 
estimate of the number of contaminants with true colour 1 .7 < 
i' —z'< 2.0 is done in two steps as explained below and is done 
separately for each patch. 

1. Six patches of the CFBDS have /-band follow-up for all 
the candidates with measured colour i! — zf > 1.7. These 
patches contain 228 candidates over a total area of 199 deg 2 . 
The zf—J diagnostic shows that 35% of these objects have 
indeed spectral type later than M8; the 65% another ob- 
jects are <M8, artefacts or quasars. We assume that this 
value of 35% is representative of the fraction of >M8 in the 
1.7 < i' - zf < 2.0 sample for all the patches. These >M8 
dwarfs are the main sources of contaminations in our >L5 
stellar count 

2. Thus we compute the probability to get a M8 to L4.5 dwarf 
with true colour 1.7 < i' - zf < 2.0 but resulting colour 
i' — zf > 2.0. We draw a photometric error in the i' and zf 
bands for 35% of the objects with 1.7 < i' —zf < 2.0. Objects 
with resulting i'-zf > 2.0 and zf < 22.5 are the contaminants. 
In practical, we added a photometric error to all objects with 
1.7 < i' - zf < 2.0, and each object that enters our selection 
box will count for only 0.35 contaminant, which provides a 
smoother contaminant distribution. 

This analysis underestimates the number of contaminants be- 
cause we do not consider direct contamination from the 1.3 < 
i' - zf < 1.7 colour range to i' - zf > 2.0, but also overestimates 
the number of contaminants because some objects within our 
1.7 <' i — zf < 2.0 contaminants sample are themselves contam- 
inants with bluer true colours. However these effects are not a 
significant problem, especially because they tend to compensate 
each other: fro m a statistical analysis of 300 000 fake objects, 
iDelormd (2008) Q derived the contamination percentage from a 
colour bin to another and showed that the overall probability for 
any of the objects in the 1.3 < i' — zf < 1-7 colour range to be 
falsely taken into account in our contamination calculation as a 
potential contaminant is 0.35%. 

The number of contaminants in our 444 deg 2 sample is ob- 
tained in zf magnitude bins of 0. 1 mag and if - zf colour bins 
of 0.1 mag. The number of objects in these small bins is usually 
smaller than 1 . This number represents the probability of con- 
tamination within one bin. The total number of contaminants is 
30, among which 23 have a resulting colour i' - zf < 2.2. 



3.4. Contaminants 

Reddened stars are not possible contaminants because the 
CFBDS fields are in a low galactic extinction area (/' - zf galac- 
tic reddening < 0.02 mag). Giant stars are not expected in our 
sample because with the apparent magnitude of our sample (z' 
between 18.5 and 22.5) they would be located at 100 kpc to 1 
Mpc, well outside of the Milky Way. Actually, we noticed a con- 
tamination by giant stars for fields close to Messier galaxies that 
we had to eliminate from the survey. 

Extremely reddened galaxies could potentially have if - zf 
and zf - J colours the same as those of brown dwarfs. These 
galaxies would have to be at a fairly high redshift (z > 1) for two 
reasons: (i) higher redshift means shorter rest-frame wavelengths 
and therefore a lower Ay to give the red observed colours; (ii) 
the CFBDS objects are unresolved in good seeing (typically 
0.6 - 0.8") and lower redshift galaxies would be easily resolved. 
However, the high redshift necessary, coupled with the fairly 
bright / magnitudes of the CFBDS objects and the fact that the 
red zf - J indicates significant extinction still at / (e.g. about 3 
magnitudes of extinction for rest-frame E(B - V) = 1), means 
these galaxies would need to have extremely high intrinsic ab- 
solute magnitudes and stellar masses (> 2 x 10 12 M o ). Galaxies 
of thi s mass are extremely rare at redshift z ~ 1 ( Drorv et alj 
I2009h . and are compo sed of evolved ste llar populations rather 
than dusty starbursts dTavlor et al.ll2009h . These massive galax- 
ies would lik ely be resolved by gr ound-based seeing out to red- 
shift z ~ 2 dMancini et alJl2010h . In conclusion, we find that 
high-redshift unresolved galaxies cannot appear as red as the L 
and T dwarfs and do not cause a significant contaminant. 

As already mentioned, high-redshift quasars have the same 
i' - zf colour as brown dwarfs but are easily removed form the 
sample through their bluer zf—J colour by about one magnitude. 
Moreover, they are rare objects, so their contamination is indeed 
negligible. 

As shown in the colour-colour graph in Fig. Q]our sample se- 
lects L5 and cooler dwarfs. However, each Megacam field shows 
on average 50000 sources, of which about one is selected as a 
brown dwarf candidate, meaning that non-brown dwarf sources 
outnumber brown dwarfs by more than a factor 10000, making 
3 or 4 <t photometric noise scattering a probable source of con- 
tamination, particularly in the i' band where brown dwarfs are 
faint and observed at low signal-to-noise ratio and a fraction of 
earlier type dwarfs with true colour i' - zf < 2.0 are sent within 
the i' — zf > 2.0 sample (see Fig. 01. This contamination has 
to be estimated carefully because earlier dwarfs are much more 
numerous in a limited-magnitude sample than the later ones. 

There is only a small probability for a given object to 
have large photometric deviations in two independent filters. 
Therefore the /-band observations at a good S/N ratio are use- 
ful to identify a significant fraction of the contaminants. All ob- 
jects with zf - J < 2.5 are removed from the sample. This selec- 
tion rejects (1) artefacts with obviously no J-band detection, (2) 
quasars, and (3) dwarfs with true i' - zf colour lower than 1.3, 
corresponding to spectral type earlier than M8, and reddening by 
photometric noise. 

The contamination in zf - J does exist but has a much lower 
impact on our study for several reasons: (i) we only obtain / 
magnitudes for the selected candidates, 249 objects, which make 
several sigma noise excursions non-significant. This is espe- 
cially true because the signal-to-noise ratio in / is much better 
than in /', leading to much smaller photometric errors, typically 
in the 0.02-0. 1 magnitude range, (ii) Because the colour spectral 
type relation is very steep at the M7-M8.5 transition, going from 



5 http://tel.archives-ouvertes.fr/docs/00/35/10/10/PDF/these.corrigee.012009. 
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4. The brown dwarf space-density 

The CFBDS patches contain a sufficient number of images and 
candidates to provide significant statistics. Completeness and 
contamination are thus handled separately in each patch before 
combining all data to derive the brown dwarf space-density. 

In order to take into account the photometric errors, each 
object is spread over a colour range following the probability 
error distribution already used for contaminants (FiglUi. In this 
process, we take into account that the true colour of the object 
is more likely to be bluer than redder, considering the density 
gradient as a function of i' - z' colour: the probability shown 
in Figj4]is weighted by the ratio of densities at the new colour 
and the observed colour, where densities as a function of colour 
are derived from our sample. Resulting objects with new i' - 
z' colours and new z' magnitudes beyond our selection limits 
are removed. Each object is thus split into several objects (about 
one hundred) whose weight is smaller than one, representing the 
actual likelihood that the object belongs to a given colour bin 
for a given observed colour. The final weight w is obtained by 
dividing by cy, the completeness in the colour-magnitude bin 
(z, J). The contaminants are assigned a negative weight. 

We determine the luminosity function of ultracool field 
dwarfs with the generalised form of the V max classical technique 
dSchmidjll968h . which allows to consider density gradients in 
the Galactic disc. We compute the maximum volume probed by 
our magnitude-limited survey at a given colour (that is at a given 
absolute magnitude). This geometric volume is corrected to take 
into account the decrease of the stel lar density with increas- 
ing distance above the G alactic plane dFeltenll 19761 [Stobie et al. 
ll989HTinnev et al.lll993h : 

V Gen = fi-r-P - {f + 2^ + 2) exp(-£)], 
sin \b\ 

where Q is the area of the survey, b is the Galactic la titude of the 
field, h is the thin disc scale height (h ~ 250pc from iRobin et all 
l200l . and £ = ds1 ^ with d the maximum distance of detection. 

For each object d is estimated. Thus the object is counted as 
the inverse of the maximum volume Veen in which it is observed. 
The luminosity function is the sum over all objects within an 
absolute magnitude bin: 

where w is the weight assigned to each object in this colour 
bin, taking into account the photometric errors and the complete- 
ness, as explained above. 

The parallax is obtained from absolute magnitude - colour 
relations as described below. 

4.1. Absolute magnitude versus colours and spectral type 
relations 

Among the publicly available spectra found in the L and T dwarf 
data archive that are used to derive synthetic colours (Sect. I3.ll ), 
32 also have measured parallax. They allow us to derive an 
absolute magnitude colour relation. The absolute magnitude- 
colour diagrams are shown in Figs. [5] and [6] The symbols in- 
dicate the different spectral classes: open circle for a M7 dwarf, 
open squares for M8 to L4.5 dwarfs, filled circles for L5 to L9.5 
dwarfs, open triangles for TO to T4.5 dwarfs and filled triangles 
for T5 and later dwarfs. The dotted line shows the derived abso- 
lute magnitude - colour relation. The scatter around the relations 
ranges from 0.2 to 0.8 mag with a mean value of 0.5 mag. 



As one can see in the upper panels of Fig. [5] the i' - z' 
colour is a good luminosity estimator for late M and L dwarfs, 
as sh own before in the SDSS photometric system dWest et alj 
2005). In the T dwarf domain, the i' - z' colour is a poor lumi- 
nosity estimator. Moreover many T dwarfs are not detected in 
the i' band. Thus we use the z' — J colour as a proxy. However 
the absolute z' and J band magnitudes of early T dwarfs are 
nearly constant. This translates into a flux increase in the z' and 
J bands from L8 to T4 spectral types (known as 7-band bump , 
iDahnet al.ll200l iTinnev & Burgasseril2003l: IVrba et alJ l2004). 
contrarily to the classical decrease of luminosity with increas- 
ing spectral type. This behaviour at the L-T transition might be 
caused by the clearing of dust clouds in the atmosphere that re- 
yeals a deeper and higher temperature re gion of the atmosphere 
( Bur gasser et al]|2002t iKnapp et al.ll2004l) . 

Of our candidates, 32 were found to be T dwarfs from their 
spectroscopic follow-up. They are plotted in Fig. [6] (circled 
symbol) with their observed z' - J colour and with My and 
Ml determined from brown dwarfs with known parallax that 
have the same spectral type. Whereas the agreement is good for 
early T dwarfs (open triangles), the synthetic colours of late T 
dwarfs may be biased and are redder than the observed colours 
of T dwarfs with similar spectral type, as already discussed in 
Sect. 13.11 We performed spectroscopic follow-up for all of our 
candidates with z' - J > 4. Thus we do not use their z' - J 
colour to derive their absolute magnitude but use their spectral 
type instead. Absolute magnitude versus spectral type relations 
are shown in Fig. [7] 

To summarise, we use 

- M z i and My vs i' - z' relations for i' - z' < 2.4 (L dwarfs): 

M z , = 3.32(i - z') + 9.86; My = 2.90(/ - z') + 7.70 

- M z > and Mj\sz' -J relations for i' -z' > 2.4 (early-T dwarfs): 

M z , = 0.74(z' - J) + 14.86; My = 13.93 

- M z i and My vs spectral type relations for dwarfs with spec- 
troscopic follow-up (part of the early-T dwarfs and all late-T 
dwarfs). 

for <T4: My = -0.28 * spT + 17.49; M z! = My + (z - J) 
for >T4: My = 0.89 * spT + 0.55; M z , = Mj + (z' - J) 

4.2. Photometric distances 

The previous relations are used to derive the photometric dis- 
tances of the brown dwarfs in our sample. The scatter in the 
absolute magnitude-colour relation translates to an uncertainty 
of 26% on the distance estimate. Figure [8] shows the distribu- 
tion of photometric distances obtained from the z' magnitudes 
without (solid line) and with (dotted line) contamination and 
completeness correction of the sample. It mostly contains brown 
dwarfs within 20 pc and 120 pc, which most probably belongs 
to the galactic disc, although it may contain close subdwarfs 
that can be revealed by high kinematics or low-metallicity fea- 
tures in the spec trum (Burgas ser et al.ll2~002tlGeballe et al.l l2002: 
Burg asser et al] l2009t iBurningham et al.ll2010al) . The distribu- 
tion peaks between 60 and 100 pc, which is about the maxi- 
mal detection distance for late-L and mid-T dwarfs. Most of the 
objects detected farther away than 120 pc are early T dwarfs 
brighter than both late L and mid-T because of the J and z'- 
bump. 
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Fig. 8. Photometric distance distribution computed in the z' 
band of our sample, without (solid line) and with (dotted line) 
contamination and completeness correction. 



4.3. Malmquist bias 

The Malmquist bias (Malmquist 1920) comes from the intrin- 
sic dispersion in the absolute magnitude-colour relations and the 
limited-magnitude definition of the sample. A given colour (or 
spectral type) does not correspond to a unique luminosity, but 
rather a distribution due to intrinsic scatter in metallicity and age, 
and non detected binaries that appear brighter for their colour. 
This dispersion has two distinct effects: 

- for a limited-magnitude sample, the mean absolute magni- 
tude observed M is lower than the true mean magnitude Mq 
of objects with a given colour: among the most distant ob- 
jects, the intrinsically most luminous ones are detected only 
(where the true intrinsic absolute magnitude Mq > M), 

- the number of stars with given apparent magnitude and 
colour is larger, the brightest stars being observed at larger 
distances d and the number of objects increasing as d 3 . 

iMalmquistl d 19201) gave the first correction method for this 
bias, assuming a Gauss ian distribution o f the luminosity at a 
given colour, and later lStobie et alJ dl989l) proposed the follow- 
ing analytic correction: 

— a (0.6/nlO)V -0.6M0cr 2 — , 
O $ 

where <x is the luminosity dispersion, <£ the luminosity function, 
and <£>' = 4%. The first term is the volume element correction, 
the second term corrects for the magnitude translation from M 
to Mq. 

Because we already took into account the uncertainties in the 
colours due to photometric errors, <x only depends on the scat- 
ter in the absolute magnitude-colour relation. This scatter ranges 
from 0.2 to 0.8 mag. We consider the mean value cr = 0.5 mag. 

4.4. The luminosity function of systems 

Our magnitude-limited sample is also biased by binarity effect. 
An object in a non-resolved binary system appears brighter and 



enters our sample contrarily to the same isolated object. This 
also affects the colour of the objec ts. Direct imaging surveys 
dGizis et al.l l2003t iBouv et all 1200 3l) that probe separations of 
>2 AU show that ~ 10 - 15% of field brown dwarfs are bina- 
ries with separations that peak between 2 and 4 AU. Lower sep- 
arations probed by spectroscopic surveys (see iJoergensI ([2008) 
and references therein) contribute a further 7^% (separations 
< 0.3 AU) or 10^ 8 % (separations < 3 AU). The true BD binary 
fraction is likely betwe en 15 and 25%, similar to estimates from 
iBasri & Reinersl d2006l) . 

The purely observational data are not the luminosity func- 
tion but a stellar count as a function of colour and magnitude. To 
obtain the luminosity function of the systems a first transforma- 
tion is done with colour-magnitude or spectral type-magnitude 
relations. Bias could be introduced during this first step, but the 
relations used (and their limitations) are relatively well known. 
So the luminosity function of the systems is a good proxy to 
compare stellar counts of different surveys (with different sets 
of filters and magnitude limits). Moreover, it is the luminosity 
function of systems which is used for estimating the number of 
sources in any survey (with different seeing and spatial resolu- 
tion). 

As explained above, the correction of this function in a true 
luminosity function (density of stars by magnitude bin) requires 
the knowledge of the statistical multiplicity (fraction of binaries, 
triple and higher order systems unresolved and distribution of 
mass ratios between the components). To date such statistics for 
brown dwarfs have still a significant margin of error. Hence in 
this paper we chose to present the luminosity function of systems 
defined as the luminosity function uncorrected from the bias due 
to multiple stars. 

Figure [9] shows the luminosity function of systems in the z' 
and J bands obtained from the brown dwarf sample corrected 
from contamination and completeness. The number of objects 
in magnitude bins are also indicated. The values are listed in 
Table [7] Error bars represent the 95% confidence level interval, 
derived from Bayesian statistics assuming a Poissonian distri- 
bution and a co njugate prior following the method described by 
iMetchev etafl d2008l) . 



Table 7. Brown dwarfs /-band and z'-band luminosity function 
<S> (10~ 3 objects pc~ 3 mag -1 ). Error bars give the 95% confidence 
level interval, derived from Bayesian statistics. 



Mj 


N M 




M-j 


N m 




13.75 


18.0 


>0.78t* 


16.75 


26.5 


0.983$ 


14.25 


52.7 


9 qc+0.60 
^•^ -0.53 


17.25 


47.7 




14.75 


21.2 


2-80!^ 


17.75 


23.3 


Q 1 Q + 1.01 

- M8 -0.83 


15.25 


8.6 


i qo+1.06 
1,7 -0.78 


18.25 


5.3 


1 75 +1 29 

L - ' -0.86 


15.75 


2.5 


T QQ + 3.63 


18.75 


0.9 


0-042J1 


16.25 


1.2 


3 75 +7 - 23 


19.25 


2.1 


0.813£ 


16.75 


0.6 


9 sq+8.44 


19.75 


1.8 


9 SQ+3.78 
iJ -1 .87 


17.25 


0.9 


2.78+f« 


20.25 


1.2 


3 16 +60S 

J - -2.55 


17.75 


1.9 


13.75!« 7 f 


20.75 


1.2 


9 17+4.I8 
A-1/ -1.76 
1 3 79+17.46 

LJ. 1 ^- g 29 








21.25 


2.2 



One notices an increase at the faintest absolute magnitudes, 
due to T8 and later dwarfs. However, the uncertainties are very 
large, based on three objects only and the luminosity function 
remains consistently flat. Furthermore, the luminosity-spectral 
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type relation used for distance estimate is rather rough (see 
Fig. |7}. Assuming a luminosity 0.5 mag brighter leads to a 
luminosity function twice smaller in these bins. However, the 
T8.5 dwarf found aroun d a M4 dwarf with measured parallax 
(Burningham et al. 2009, open square in Fig. |7]i tends to show 
that the late-T dwarfs luminosity is not underestimated, on the 
contrary. This uptick in numb er count is expe cted from luminos- 
ity function simulations (see lAllen et al.l2005l and the pile-up of 
objects with T e ff < 500 K in their Fig. 2). 

Finding half of the known ultracool brown dwarfs (>T8) 
could have been a statistical fluke. Future ultracool brown dwarf 
discoveries will build up this statistic and check whether this sig- 
nificant increase at the lower end of the luminosity function that 
we find here is real or not. If true, this would mean either that 
the lowest mass brown dwarfs are more numerous than expected 
- which would be different t o what is observed in young cluster 
(see e.g. lMoraux et al. I l2007h - or that the number of old brown 
dwarfs of any mass is high. Indeed, whatever their mass, brown 
dwarfs cool down to the latest and coldest end of the spectral 
range after a few Gyr and should accumulate there because the 
cooling rate significantly slows down at these low temperatures. 
This speculation w ould agree interestingly with the hints from 
iGould et"al] ((2009) that old-population brown dwarfs could be 
much more numerous than young ones, based on their a priori 
low probability detection of a thick-disk brown dwarf in a mi- 
crolensing event. 

5. Discussion 

Table [8] summarises the comparison with other br own dwarf 
space densities in My bins: the space density from ICruz et al.l 

(2007) computed from a 20 pc volume-limited sam ple of M7 to 
L8 dwarfs and the one obtained I Allen et al] d2005l) from a sam- 
ple of 14 T-dwarfs and a volume-limited sample of late M and L 
dwarfs. As seen in Fig.[TJ the T dwarfs can be selected from their 
i' -z' colours larger than 2.6. For objects with no i' detection, the 
z' — J > 3.5 limit is also a good selection criterion for T dwarfs. 
Thus it is possible to compute the space density for T-dwarfs in 
differe nt spectral type ranges and to compare with lMetchev et all 

(2008) result . Our results agree with other s. For early-T - dwarf s, 
they agree with the lower value found bv lMetchev et al.l ([2008). 

The most important aim in deriving the luminosity function 
for ultracool dwarfs is setting constraints on the mass function. 
The computation of mass function from the luminosity function 
is quite straightforward for stars whose luminosity follows a sim- 
ple function of their mass. The case is more complex for brown 
dwarfs that undergo continuous cooling and gravitational con- 
traction during their life. The decrease of luminosity with time 
provokes a degeneracy: a low-mass brown dwarf can be as bright 
as a higher mass brown dwarf if younger. The only way to de- 
termine age and mass is to compare the spectra with theoretical 
atmosphere models, which still give quite uncertain results, be- 
cause the physics of these cool atmospheres is very complex and 
not yet totally understood. The non-relation between mass and 
luminosity makes difficult the determination of the mass func- 
tion. A detailed determination of the substellar mass function is 
beyond the scop of this paper. 

However, first estimates of the subst ellar mass functio n have 
been obtained in the L-dwarf domain: iReid et all d 19991) with 
1 < a < 2, lBurgasserid2002l) 0.5 < a < l. lAllen et alJd2005l) a < 
0.3 + 0.6, where a is the slope of the power-law mass function 
T'(M) = dN/dM oc M~ a . These results are consistent with the 
subst ellar mass function in young open clusters showing a 0. 6 
(e.g. ICaballero et all 120071: Moraux et al.ll2007b lLuhmadl2007h . 



Recently, Met chev et alJ d2008l) derived a T dwarf space density 
that was mostly consistent with a — an d the compar i son of 
the observed number of T4-T8.5 dwarfs by Pin field et alJ (2008) 
favoured a < 0. The analysis of 47 T-dwarfs found in the Large 
Area Survey (LAS) of UKIDSS also suggests that the substellar 
mass function is declining at lower masses ( Burningham et al. 
I2010bl) 

iBurgasserl J2007) performed simulations assuming different 
mass functions. The comparison between our luminosity func- 
tion and these simulations is shown in Table [9] It suggests that 
the best agreement is obtained with a flat mass function (a = 0) 
or even a decreasing mass function. Note that our values dis- 
agree for early T-dwarfs. However, the physics at the L/T tran- 
sition is difficult to model, not only due to the cooling of the 
atmosphere but also to the clearing of the atmosphere that needs 
to model properly th e hydrodynamics and the clouds formation 
dFrevtag et al.l 12010). making the t imescales uncer t ain. If con- 
firmed, this low value would back Burgasseret al. (2007) sug- 
ges tion that the L/T tran sition occurs rapidly. 

iPinfield etalJd2008l) suggested that a single power-law expo- 
nent is not optimal when describing the substellar mass function 
in the field, or that the different measured a values result from 
different models used to convert between mass and magnitude. 
Further studies based on more extended samples, such as the one 
we will be able to define at the completion of the CFBDS survey, 
or from the ongoing Large Area Survey performed by UKIRT, 
are needed to make reliable investigations on the mass function 
of field brown dwarfs. 

Table 9. Comparison of t he brown dwarf sp ace densities p (10 3 
objects pc~ 3 ) obtained by Burgasser (2007) from simulations as- 
suming different mass functions: *F(M) = dN/dM oc M~ a . 



Spectral type 


p(a = 0) 


p(a = 0.5) 


p(a = 1) 


p(CFBDS) 


L5-L9.5 


2.3 


2.9 


4.9 


2 +u i! 


T0-T5.5 


3.3 


4.4 


6.1 


1 4 +03 


T6-T8 


6.4 


9.5 


14.4 





6. Conclusions 

At mid-course of the CFBDS, we define a uniform sample to 
investigate the field brown dwarf luminosity-function. We ob- 
tained spectroscopic follow-up for most of T-dwarf candidates. 
Because it is not realistic to obtain spectra for all L-dwarf candi- 
dates, we used the noise properties of the images to compute 
the contamination and completeness with good enough accu- 
racy to properly characterise the sample. The sample covers a 
large spectral type range, including the L-T transition. It con- 
tains 56 L-dwarfs cooler than L5 and 41 dwarfs along the whole 
T-dwarfs sequence. We have measured the density of late L5 to 
TO dwarfs to be 2.0+°® x 10" 3 objects pc 3 , the density of T0.5 
to T5.5 dwarfs to be 1-4+q 3 , x 10~ 3 objects pc~ 3 , and the density 
of T6 to T8 dwarfs to be 5.3+I* x 10 3 objects pc 3 . Three lat- 
est dwarfs at the boundary between T and Y dwarfs give the high 
density 8.3+?jX 10~ 3 objects pc~ 3 , although the uncertainties are 
very large and the luminosity function is still consistent with be- 
ing flat. Since then, at least three ultra cool brown dwarfs (>T8) 
have also been discovered in UKIDSS (Burningham et al. 2008, 
120091) that currently probes a slightly higher but comparable vol- 
ume for late-T dwarfs. Even if these combined statistics still deal 
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with a low number of objects, this could indicate that the number 
of ultracool brown dwarfs found by the CFBDS is not entirely 
due to a fluke of statistics and that many of these objects re- 
main to be discovered in the solar neighbourhood. This can be 
expected as, whatever their mass, brown dwarfs cool down and 
finally accumulate in this spectral type range. 
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Table 5. Photometry of L5 and later dwarf candidates in the MegaCam photometric system. For objects without detection in the 
/'-band, a lower limit is given for the i' - z' colour. The spectral type is given for dwarfs with spectroscopic follow-up. 
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Warren et al. (2007). b Delorme et al. (2008a). c Chiu et al. (2006). 
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1 Chiu et al. (2006). Knapp et al. (2004). 



Table 8. Com parison of the brown dwarf space de n sities p (10 3 objects pc 3 ) obtained from CFBDS given separately for L and 
T-dwarfs, and All en et al.1 (l2005h : ICruz eTafl (|2007); Metchev et al. (2008|). 



Mi p(Cruz) p(Allen) p(Metchev) p L (CFBDS) p r (CFBDS) p M ,(CFBDS) 



13.75 


0.5 ± 0.2 


1.0 ± 


0.3 






0.78*$ 


14.25 


> 0.7 ± 0.2 fl 


> 1.7 i 


: 1.4° 




1 09 +0 30 
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> 0.3 ± 0.2° 


>2.3d 


: 1.4° 




9 i q+0.78 


0.62^ 
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> 1.9 i 


:0.9 fl 




fl 7Q+0.49 
' 7 -0.34 


1 14 +112 

1-1 -0.67 
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2.0 ± 
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57 +1 63 

"■-"-0.52 


9 <T9+3.67 
A - J -1.82 


16.25 




4.7 ± 


0.3 






3 75 +7 - 23 


T0-T5.5 








9 Q+0.9 
^■ J -0.9 




1 4+0.3 


T6-T8 








4 7+ 3.i 

^•'-2.8 






T8.5-T/Y 












8 - 3 -?i 



" a lower limit is given when the authors noted the incompleteness of their sample due to colour selection biases. 
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Fig. 5. Mj and M' z absolute magnitudes versus synthetic colour i' - z' from available spectra of brown dwarfs with known parallax. 
Open circles: dwarfs earlier than M8. Open squares: M8 to L4.5 dwarfs. Filled circles: L5 to L9.5 dwarfs. Open triangles: early T 
dwarfs. Filled triangles: late T dwarfs. The dotted line shows the derived absolute magnitude colour relation valid for late-M to L 
dwarfs. Error bars are indicated when larger than the symbol. 
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Fig. 6. Same as Fig.|5]for the z' - J colour. Additional objects are shown (circled symbols). They are our candidates with spectro- 
scopic follow-up. For these objects, z' - /is the observed colour and the absolute magnitude is derived from their spectral type. The 
dotted line shows the derived absolute magnitude-colour relation valid for the early-T dwarfs. Error bars are indicated when larger 
than the symbol. Note that only objects with measured parallax - non circled objects - are used to derive the relation. 
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Fig. 7. Mj and M' absolute magnitudes versus spectral type from available spectra of brown dwarfs with known parallax. The dotted 
line shows the derived absolute magnitude spectral type relation for the early-T and late-T dwarfs. Error bars are indicated when 
larger than the symbol. The open square with error bars shows the T8.5 companion to Wolf 940A, a M4 dwarf (Burningham et al. 
2009). 
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Fig. 9. Luminosity function of systems computed in the z' band (left) and the J band (right) of our sample. 



